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Abstract 

The direct observation of gravitational waves (GW) in the near future, and 
the corresponding determination of the number of independent polarizations, 
is a powerful tool to test general relativity and alternative theories of gravity. 
In the present work we use the Newman-Penrose formalism to characterize 
GWs in quadratic gravity and in a particular class of f{R) Lagrangians. We 
find that both quadratic gravity and the f{R) theory belong to the most 
general invariant class of GWs, i.e., they can present up to six independent 
polarizations of GWs. For a particular combination of the parameters, 
we find that quadratic gravity can present up to five polarizations states. 
On the other hand, if we use the Palatini approach for f{R) theories, 
GWs present only the usual two transverse-traceless polarizations such as 
in general relativity. Thus, we conclude that the observation of GWs can 
strongly constrain the suitable formalism for these theories. 
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1. Introduction 

Modifications of Einstein's gravity have been considered in several 
approaches. The Lagrangians which consider higher orders of the Ricci 
scalar R and the Ricci tensor R^i, have been proposed as extensions of 
general relativity [Ij. The semiclassical theory considers the backreaction of 
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quantum fields in a classical geometric background [21 [3] . Such Lagrangians 
predict field equations with four orders derivatives of the metric, rather than 
the two orders derivatives in the general relativity theory |2l H] . Quadratic 
Lagrangians have been used also to yield renormalizable theories of gravity 
coupled to matter [S]. Moreover, higher-derivative theories arise as a low 
energy limit of string theories [HI [T| . 

Starobinsky j;8j argues that the higher order terms could mimic a 
cosmological constant. Recently, this idea has been largely studied as a 
potential way to address the dark energy problem. In this context, several 
different forms of the modified Lagrangians have been recently considered 
constituting a class of theories, the so-called f{R) theories (see e.g. and 
references therein). 

In the context of gravitational waves (GWs), it was shown that quadratic 
gravit}|^ presents a frequency-dependent shift in the wave amplitude with 
respect to the wave amplitude obtained from general relativity [10]. For the 
f{R) theories, it was found that GWs can have a massive-like scalar mode 
besides the usual transverse-traceless modes [TT] . 

A powerful tool to study the properties of GWs in any metric theory of 
gravity was developed by Eardley et al. [12] . The basic idea is to analyze all 
the physically relevant components of the Riemann tensor R\^ku-, which cause 
relative acceleration between test particles. The GWs in a metric theory 
involves the metric field g^y and any auxiliary gravitational fields that could 
exist. But the resultant Riemann tensor is the only measurable field. 

In their work, Eardley et al. used a null-tetrad basis in order to calculate 
the Newman-Penrose [13] quantities in terms of the irreducible parts of Rx^nu, 
namely, the Weyl tensor, the tracelless Ricci tensor and the Ricci scalar. 
This analysis showed that there are six possible modes of polarization of 
GWs in the most general case, which can be completely resolved by feasible 
experiments. Thus, it is possible to classify a given theory by the non-null 
Newman-Penrose quantities [T2] . 

The aim of the present work is to characterize GWs for a particular 
class of f{R) gravity and for quadratic gravity making use of the Newman- 
Penrose formalism. Since the field equations derived from such Lagrangians 
in the metric formalism yield dynamical equations for R and Rfj,u, the analysis 



^For the sake of definiteness we will call quadratic gravity that theory which takes into 
account not only i?^ but also R^j^R'^'' in the Lagrangian. 
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consists basically to find the resultant expressions for R and R^^i, in the week 
field limit, without the need to write explicitly the expressions in terms of 
the metric perturbations. This makes the classification easier and clear. 
We also mention how the use of the Palatini approach in deriving the field 
equations for f{R) gravity affect the number of independent polarizations 
of GWs. This comes from the fact that the classification of a given f{R) is 
formalism-dependent . 

It is argued that the observations of the GWs in the near future (for 
the current status of GWs detectors see e.g. [HI [151 [HI [13 [H]), and 
the corresponding determination of all possible states of polarization, is a 
very powerful test for the present studied alternative theories of gravity. 
Particularly, we show that GWs experiments can be decisive for quadratic 
gravity and in the determination of the suitable formalism for f{R) theories, 
i.e., the use of the metric or the Palatini approaches, since the number of 
polarizations of GWs depends on the formalism usecf] 



2. The Newman-Penrose formalism - an overview 

Throughout this paper we consider GWs propagating in the +z direction. 
So, all the quantities are functions only of t and z. 

At any point P, the null complex tetrad (k, 1, m, m) is related to the 
Cartesian tetrad (e^, e^;, e^^, e^) by: 

k=^(et + e,) (1) 

l=^(e*-e.) (2) 
m=^(e^ + iej,) (3) 

^=^iGx-iey) (4) 
It is easy to verify that the tetrad vectors obey the relations: 

— kl = mm=l (5) 



^See, in particular, |19j for an application of the Ncwman-Pcnrosc formalism to 
determine the GW polarizations in massive gravity. 
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km = km = lm = lm = 



(6) 



The null-tetrad components of a tensor T are written according to the 
notation: 



where (a, 6, c, ■ ■ ■) run over (k, 1, m, m) and (/i, z/, A, ■ ■ ■) run over {t, x, y, z) 
since we are working in cartesian coordinates. 

In general the Newman-Penrose quantities, namely the ten \E''s, nine $'s, 
and A, which represent the irreducible parts of the Riemann tensor Rx^nf, 
are all algebraically independent. When we restrict ourselves to nearly plane 
waves, however, we find that the differential and symmetry properties of 
Rx^iKu reduce the number of independent, nonvanishing components, to six. 
Thus, we shall choose the set {\l/2, ^E's, ^^4, $22} to describe, in a given null 
frame, the six independent components of a wave in the generic metric theory. 
In the tetrad basis, the Newman-Penrose quantities of the Riemann tensor 
are, therefore, given by: 



Note that, \E'3 and \E'4 are complex, thus each one represents two independent 
polarizations. One polarization for the real part and one for the imaginary 
part, totalizing six components (see fig. [l] which was taken from [20]). 

Analyzing the behaviour of the set {\E'2, ^3, ^4, $22} under rotations, we 
see that they have the respective helicity values s = {0, ±1, ±2, 0}. 

We also have the very useful relations for the Ricci tensor: 




(7) 



^2 — —-p^Rlklk, 


^3 = —^Rlklm, 
^4 = —Rmm, 

$22 = —Rlmlfn- 



(8) 



(9) 
(10) 

(11) 



Rlk — Rlklk, 



(12) 



Rll — '^Rlmlm, 

Rim Rlklrrii 
RlWi Rlklrni 



(13) 
(14) 
(15) 



and for the Ricci scalar: 



R — —'^Rik — —'^Rlklk- 



(16) 
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Figure 1: The six polarization modes of weak, plane, null GW permitted in any metric 
theory of gravity. Also shown is the displacement that each mode induces on a sphere 
of test particles. The wave propagates out of the plane in (a), (6) and (c), and it 
propagates in the plane in (d), (e) and (/). The displacement induced on the sphere of 
test particles corresponds to the following NewmanPenrose quantities: Re^4(a), Im^4(6), 
$22(c), *2(d), Re*3(e), Ini*3(/). 



The six amplitudes {\l/2, ^l/s, ^1/4, $22} of a wave are generally observer- 
dependent (for more details see the Apendix |A]). However, there are 
certain invariant statements about them that are true for all standard 
observers if they are true for any one. These statements characterize invariant 
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£■(2) classes of waves. The name of each class is composed of the Petrov type 
of its nonvanishing Weyl tensor and the maximum number of nonvanishing 
amplitudes {\E'2, ^E's, ^1/4, $22} (the dimension of representation) as seen by 
any observer. Both the Petrov type and the dimension of representation 
are independent of observer. Considering standard observers such that: (a) 
each observer sees the wave travelling in the -\-z direction, and (b) each 
observer measures the same frequency for a monochromatic wave, then the 
£^(2) classes in order of decreasing generality are: 

• Class 11^: \E'2 7^ 0. All standard observers measure the same non-zero 
amplitude in the ^2 mode. But the presence or absence of all other 
modes is observer- dependent; 

• Class Ilh: \&2 = 0, ^37^ 0. All standard observers measure the 
absence of \E'2 and the presence of But the presence or absence of 
^'4 and $22 is observer- dependent; 

• Class A^s: vl/2 = vE^g = 0, ^4 ^ 0, $22 7^ 0. Presence or absence of all 
modes is observer- independent; 

• Class A^2: *2 = *3 = ^22 = 0; ^4 7^ 0. Observer-independent; 

• Class Oi. \Ef2 = ^3 = ^4 = 0; $22 7^ 0. Observer-independent; 

• Class Oq: ^2 = ^3 = ^4 = $22 = 0. Observer-independent. All 
standard observers measure no wave. 

3. Polarization modes of gravitational waves in f{R) theories 

3.1. The metric formalism 

Let the gravitational action be an arbitrary function of the Ricci scalar: 



By varying this action with respect to the metric g^^i, we have the following 
vacuum field equations: 




(17) 



(18) 



where the prime represents derivatives with respect to R. 
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In the sequence we restrict ourselves to the following class of the f{R) 
theories: 

f{R) = R- aR-^ (19) 



Substituting (19) in the field equations (18) we obtain the following relations 



between R^^ and the Ricci scalar R: 



{R - aR-P)g,, + 2al3V ,R-^^+^^ - 2a(3g,,UR-^^+P) 



2[1 + 

Contracting this expression we have a dynamical equation for R: 

R = a[{f3 + 2)R-'^ + 3/?ni?-(^+^)] 



(21) 



The classification procedure involves examining the far-field, linearized, 
vacuum field equations of a theory. In what follows we examine different 
cases for the possible values of a and (3. In each case, we first find R from 
equation (21) and then we can compute R^u from (20). 



• Case a = 

This is the trivial case, which reduces to the general relativity theory. 
From the equations (21) and (20) we find i? = and R^^ = 0. Consequently, 
from the relations of the section [2] we deduce that: 



R, 



Iklk 



R, 



imlm 



R 



iklm 



R, 



Iklm 







and so we have: 



^2 = ^3 = '^22 = 0. 

And since we have no further constrains: 

^4 ^ 0. 

And, as expected, the E{2) classification for general relativity is N2. 



(22) 
(23) 

(24) 



• Case a^0,(3>l 
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For a 7^ 0, the equation (21) can be written as: 



□i?-(i+/3) + ^R-^ - -1-R = 0. 



(25) 



3p 3ap 

Working in the weak field regime, if /3 > 1 we have R^^ ^ R and the 



equation (25) now reads: 



+ 2 P 
□0+^0^ =0, 

where we have renamed = R^^^^^\ But this equation is of the form: 

U(h- — =Q. 



(26) 



(27) 



Comparing (26) and (27) we have the potential: 

•(/3 + 2)(/? + l)" 



3/5(2/5 + 1) 



2/3+1 
) (3+1 



(28) 



Since the field is Lorentz- invariant, we can solve the equation (27) by a 



(27), i.e., the solution of the equation: 



very known method (see, e.g., [2T]). First, consider the static solution of 

(29) 



d?(t) _ dU 
dz'^ d(j) ' 



which can be written as: 



1 

2 V dz 



(30) 



Assuming the potential (28) in (30) and integrating we have: 

2(/3+l) 



where 



2(/3 + 2)(/3 + l) 



2(/? + l) 



3/?(2/? + 1) 



1/2 



(31) 



(32) 
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and 00 = (Pi^o) is the value of the field in some initial position zq. Now, 
because the system is Lorentz-invariant, given the static solution (31), one 
can Lorentz-transform it to obtain the time dependent solution: 



And, therefore, the Ricci scalar reads: 

[z - Zo) - vt 



2{/3+l) 



(33) 



R{z,t) 



1/2 



(34) 



where v is the wave propagation velocity. This is the solution of the equation 



(26) as can be verified by direct substitution. 



Substituting the solution (34) in (20) we find the non-zero components of 

2{/3 + 2y 



the Ricci tensor to first order in R: 



Rtt 



6p 



(1 - 2/?) 



Rtz — '-^R 



R, 



1 

6/3 



2(/3 + 2) 



{2(3 - 1) 



Therefore, from the relations of the section 121 we find 



R. 



R 



iklm 



R, 



Iklm 



0, Rlklk 7^ 0, Rlmlm 7^ 0, 



and so: 



^2 ^0, ^3 = 0, ^4^0 and $22 ^ 0, 
and the E{2) classification for this case is IIq. 

• Case a ^ 0,/3 < -2 



(35) 

(36) 
(37) 

(38) 
(39) 



Considering now /? < — 2 we have R ^ <^ R and the equation (25) reads: 

1 



□c 



3a/5 



-l/(l+/3) ^ 0^ 



(40) 
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where has the same definition presented above. However, now the potential 
is given by: 

Sap 

from which we obtain the static solution: 



/3+2 
2(/3+l) 



13+2 



where: 



c 



P + 2 



{(3 + l)^/MP' 

And after a Lorentz transformation we find the full solution: 



[z,t) 



' (z-Zo)-vt 5^ 



^ 2(/3+l) 
/3+2 



Thus, the evolution of the Ricci scalar for this case is: 



(42) 



(43) 



(44) 



R{z,t) 



■ (z- zo)-vt ^^-i±l' 



(3+2 



The components of the Ricci tensor reads: 



Rtt 



R 



{(3 + l){l-v^) 2 

(3 V 



R. 



tz 



R. 



1 



3(/? + l)(l-t;2) 



[3 



R, 



1 1 

+ 



R. 



(45) 



(46) 

(47) 
(48) 



3 [{f3 + 1) {1 ~ v^) 2_ 
while all the others components are null. 

Therefore, together with the formalism presented in the section [2| we can 
deduce that: 

^2 ^0, ^3 = 0, ^4^0 and $22 ^ 0. (49) 
And we are lead to classify the theories with /5 < — 2 in the E{2) class J/g. 
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• Case a ^0,(3 = -2 

This is a particular case, where the behaviour of the Ricci scalar and Ricci 



tensor are oscillatory. That is, if we use P = 
obtain: 

OR R = 0, 

6a 



-2 in the equation (25), we 

(50) 



with the solution: 



R = Roexp^ikax"") 



1 

6a 



(51) 



Considering this solution in the equation (20) with i? ^ 1 we find the non- 
null components of the Ricci tensor: 



R^^ = -(4afc2 - 1)R 



Rtz 
R. 



-2ak\lk^ 
1 



6a 



-R 



6 



(12ar + 1)R 



So, again we have: 



^2 ^0, ^3 = 0, ^4^0 and $22 ^ 0, 



(52) 

(53) 
(54) 

(55) 



and the E{2) classification is IIq. 

As can be seen, for all the studied cases (except the case a = ), the 
theory given by eq.(19) is classified in the class JJg, i.e., the most general 
classification where all the six polarization can appear for some specific 
Lorentz observer, but the amplitude is observer-independent. 

3.2. The Palatini approach 

In the Palatini approach, the metric g and the (usually torsionless) 
connection F are considered as independent variables, entering the definition 
of the Ricci tensor. The vacuum field equations, derived from the Palatini 
variational principle applied in the action (17) are: 

1 



f'R 



(p.v) 



-J 9, 



0, 



(56) 
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^liV^f'gn = 0, (57) 

where Vj^ is the covariant derivative with respect to T. We shall use the 
standard notation denoting by R(^u) as the symmetric part of R^^, i.e. 
R(^fj,i,) = \ {Rfiu + Rv/i)- It "was shown that the vacuum field equations (56) 



leads to 'universal' equations for a wide range of functions f{R) [22j. These 
universal equations are just Einstein equations with cosmological constant 
A. 

Thus, the properties of vacuum GWs in the f{R) gravity using the 
Palatini approach reduces to the problem of GWs in the Einstein equations 
considering the cosmological term. In a recent work, Naf et.al |23j| have 
analyzed this case. They expanded the perturbations in a de Sitter and an 
anti-de Sitter background. Since the Minkowski metric is not a solution of the 
vacuum field equations this approach seems to be the most straightforward. 
Considering terms up to linear order in A they calculated the non-null 
components of the Riemann tensor and found that A does not introduce 
additional polarization states for the GWs. Moreover, they shown that the 
cosmological term introduces tiny modifications in the amplitude of the wave 
which are well below the detectability of the present GWs detectors. 
Therefore, we can conclude that GWs in the f{R) gravity using the Palatini 
approach has only the two usual polarizations of general relativity, i.e., the 
polatization + and x. 

4. Polarization modes of gravitational waves in quadratic gravity 

For completeness we analyze the polarizations for the following 
Lagrangian of quadratic gravity: 

1 = j rf^xv^ [R + aR^ + iRf.uR^''] . (58) 

Variation of this action in respect to the metric gives the field equations: 

G^^ + aH^^ + 7V = 0, (59) 



where: 



H^,, = 2V^V,i? - 2g^,UR + ]^g^,R^ - 2RR^, (60) 



and 



lfj.u = y fNuR — -^g^iuOR — D-R^j. — 2R^°'Rau + ■^gfj.uRaisR"'^ (61) 
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The field equations (59) were analyzed in the linearized regime by de Rey 
Neto et. al [in] using a perturbative approach developed from the concept 
of regular reduction of a system of partial differential equations. Working 
only with the transverse-traceless part of the metric perturbations hj^ , they 
found a frequency- dependent correction in the gravitational wave amplitude 
due to the presence of the Ricci-squared term in the gravitational action. 
Here, we do not write explicitly the field equations in terms of the 
perturbations. Instead, we consider the Ricci tensor and the Ricci scalar 
as first order functions of the general metric perturbations h^p. Then, we 
find solutions for the dynamical equations to linear order in R^y and R. 
These solutions enable us to find the non-null Newman-Penrose quantities 
and classify the quadratic gravity in analogy with was made for f{R) theories 
in the preceding section. 

Considering the equation (59) to linear order in R and R^u we have the 



following equation for the Ricci tensor: 

1 



-R 



7 



-S, 



7 



where: 



{2a + 7) 



df,dyR 



Vf^uR 



4(3a + 7) 



Taking the trace of (62) we find: 



2(3a + 7)ni? + i? = 0. 



(62) 



(63) 



(64) 



From this equation, we can find that a particular combination of parameters, 
namely 7 = —3a, leads to i? = 0. For this case the solution of equation (62) 
reads: 

1 



R 



exp(igax"j 



7 



Since there are no further constrains we find: 



Rlklk — 0, Rlmlm 7^ 0, Rlklm 7^ 0, Rlklm 7^ 0. 

Therefore: 

^2 = 0, ^3 7^ 0, ^4^0 and $22 ^ 0, 

and the correspond E{2) class is III^. 



For the case 7 7^ —3a, the equation (64) gives: 



R = Roexpiiklx''), klk"^ 



2(3a + 7)' 



(65) 

(66) 
(67) 

(68) 
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With this expression in (62), the full solution for the Ricci tensor reads: 
where: 



A'^^ = -i?o(3a + 7) 



^1 ^1 ~^ 



n 



4(3a + 7) 



and: 



2(3a + 7)' 



7 



Now, from (69) we write explicitly all the components of -R^j^: 



Rtt = I [4(3a + 7)^2 - 1] i^oe'^'^i'""*) + Bue'^''''-^'^ + c.c, 
6 



(69) 
(70) 

(71) 
(72) 



Ru = -^(3a + 7)^X' + 2(^ \ 3^) -Roe^^^"^~"*^ +i?t.e-^"^^-"*^ +C.C., (73) 



1 + ^(3a + 7)w^ 



and all the other components satisfies: 



(75) 



where j = x, 2;. 

Therefore, since there are no further constrains on B^^, all the components 
of the Riemann tensor in the tetrad basis are non-null: 

Rlklk 7^ 0, Rlmlm. 7^ 0, Rlklm 7^ 0, Rlklni 7^ 0, (76) 

and so, all the Newman- Penrose quantities are also non-null: 

^2 ^ 0, ^3 ^ 0, ^4 7^ and $22 ^ 0. (77) 

Thus, the E{2) classification for the quadratic gravity in the most general 
case is IIq. 
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5. Conclusions 

Recently, working in the metric formalism, Cappoziello et al. [H] showed 
that GWs in f{R) gravity can have a massive-like scalar mode and a 
longitudinal force besides the two polarizations which appear in general 
relativity. This agrees with our result that the quantities $22 and are non- 
null for the particular class f{R) = R + aR~^ . The $22 and \E'2 amplitudes 
correspond, respectively, to a perpendicular scalar mode (breathing mode) 
and to a longitudinal scalar mode. However, it is worth emphasizing that 
since the \E'2 mode are non-null, the E{2) classification is JJg (see section 
|2|. So, for this class, it is always possible to find a Lorentz observer who 
measures all the six polarization states. On the other hand, GWs in the 
Palatini approach have only the two usual polarizations states such as general 
relativity. 

Furthermore, the method we use is not only simple but very robust and we are 
able to obtain some important information regarding the theories considered 
here. The key observational GW amplitude is \E'2. If the amplitude \E'2 would 
be detected, the f{R) models in the metric formalism like the one considered 
here would be supported and so the quadratic gravity for 7 7^ —3a. In 
this case we need another method to distinguish the two theories, we could 
compare the wave form in both cases, for example. If the mode would be 
detected, but not the amplitude \Ef2, only the quadratic gravity for 7 = —Set 
would be supported. 

Therefore, if we would be able to detect GWs, an important way to identify 
the theory of gravity could be estabhshed In the particular case of f{R) 
gravity in the Palatini approach, we showed that the polarizations of GWs 
are the same of general relativity. However, it is worth stressing that other 
information contained in the GW signals, like waveform and phase of the 
signal, could be important to permit, together with the polarizations, a clear 
identification of the theory. 
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A. Behavior of the NP amplitudes under Lorentz transformations 

We can understand the E{2) classification scheme by analyzing the 
behavior of the Newman-Penrose (NP) amplitudes {\l/2, \l/3, \l/4, $22} under 
a Lorentz transformation of the complex tetrad basis. This point is well 
explained in the reference [12]. Here we summarize the main idea. 
Consider two standard observers O and O' with tetrads (k, 1, m, n) and 
(k', 1', m', n') respectively. If we choose the k as proportional to the wave 
vector (following the convention of Eardley et al.), so we have k = k'. The 
most general proper Lorentz transformation relating the tetrads that keep k 
fixed is: 

k' = k, (78) 

m' = e^'^(m + ak), (79) 

m' = e-*'^(m + ak), (80) 

1' = 1 + am + am + aak, (81) 

where a is an arbitrary complex number which produces null rotations 
(particular combinations of boosts and rotations), while (p, which runs from 
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to 27T, is an arbitrary real phase that produces a rotation about e^. 
The transformations induced on the amphtudes of a wave by {ip, a) is: 



^2 = ^2, (82) 

^'3 = e-*'^(^3 + 3a^2), (83) 

^; = e-2^^(^4 + 4a^3 + 6cf2^2), (84) 

$'22 = $22 + 2cr\I>3 + 2a'^3 + 6o-(f\[^2- (85) 

Now, it is evident from this set of equations that the amphtudes 
{\l/2, ^1/3, ^1/4, $22} cannot be specified in an observer-independent manner. 
For example, suppose that the observer in O measure a wave having the 
only nonvanishing amplitude \l/2 [s = 0). The observer in O', in relative 
motion with respect to O, will conclude that the wave has the nonvanishing 
amplitudes \l/2, "^3, and $22 (s = 0, ±1, ±2,0). However, there is a set of 
invariant statements which define the E{2) classification scheme. Thus, we 
classify waves in an E{2) invariant manner by uncovering all representations 



of E{2) embodied in equations (82 ) — (85 ). Each such representation, in which 



some of the NP amplitudes vanish identically, is a distinct invariant class. 
The description of each class can be found in the main text. 
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